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Many important sets of normalized states in a multipartite quantum system of finite dimension 
d, such as the set S of all separable states, are real semialgebraic sets. We compute dimensions of 
many such sets in several low-dimensional systems. By using dimension arguments, we show that 
there exist separable states which are not convex combinations of d or less pure product states. For 
CO ' instance, such states exist in bipartite M®N systems when (M — 2) (TV — 2) > 1. This solves an open 

problem proposed by DiVincenzo, Terhal and Thapliyal about 12 years ago. We prove that there 
exist a separable state p and a pure product state, whose mixture has smaller length than that of p. 
We show that any real p G S, which is invariant under all partial transpose operations, is a convex 
sum of real pure product states. In the case of the 2 ® TV system, the number r of product states 
, can be taken to be r — rank p. We also show that the general multipartite separability problem 

can be reduced to the case of real states. Regarding the separability problem, we propose two 
conjectures describing 5 as a semialgebraic set, which may eventually lead to an analytic solution 
in some low-dimensional systems such as 2 g) 4, 3 ® 3 and 2 <g> 2 ® 2. 

PACS numbers: 03.65.Ud, 03.67.Mn, 03.67.-a 
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£ I. INTRODUCTION 

Q-i Entanglement reveals a fundamental difference between the quantum and classical world, which may be detected e.g. 
by Bell inequalities Q. It plays an essential role in quantum information processing, such as quantum teleportation, 
computing and cryptography. It is a hard problem to decide whether a given quantum state is entangled. We shall 
propose a new method of atacking this problem in low-dimensional quantum systems, see Conjectures [4] and [5j 

We consider a finite-dimensional multipartite quantum system described by the complex Hilbert space H = Hi <2> 
■ • • £3 Tin- We set d = d\d% ■ ■ ■ d„ where d{ = Dim Hi. For brevity, we may refer to this system also as an d\ ® ■ ■ ■ ® d n . 
' A (normalized) state is a positive semidefinite linear operator p : % — > % with Trp = 1. If the condition Trp = 1 
CO . is omitted and p ^ 0, we shall say that p is a non-normalized state. For convenience, we prefer to work with non- 
\q ' normalized states. We shall mention explicitly when we require the states to be normalized. A pure product state is 
the tensor product <£>•••<£> IV'nXV'ral) where \tpi) € Hi are nonzero vectors. A state is separable if it is a finite 

CNJ ' sum of pure product states. A state is entangled if it is not separable. There exists a necessary and sufficient condition 



> 
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for separability |ll| which shows that the separability problem is equivalent to another hard problem about positive 
linear maps. It has been shown that, in general, the problem of deciding whether a state is entangled is NP hard 0. 

Let D be the set of all normalized states on T-L, and S the subset of all normalized separable states. Let H denote 
the real vector space consisting of all Hermitian operators on W. Denote by H' the affine subspace of H defined by the 
equation Trp = 1. Our starting point is the well known fact that V is a real semialgebraic subset of H' , see Section 
ILT1 From now on we shall say "semialgebraic" instead of "real semialgebraic" . For the definition and examples of real 
algebraic and semialgebraic subsets of R" see e.g. [!, Section 2.1]. Let us just mention that H r := {p 6 H : rankp < r} 
is a real algebraic set, and so is the set H' r = H' C\H r . Since finite intersections of semialgebraic sets are semialgebraic, 
it follows that the sets T> r = T> n H r are semialgebraic. We define the length, L(p), of any separable state p to be the 
least integer r such that p can be written as a sum of r pure product states. In Proposition^ we prove that the sets 
S' r := {p E S : L(p) < r} are semialgebraic. From this result we deduce (see Corollary [3]) that S is a semialgebraic 
set, and so are the sets S r = S f] H r . 

Any closed semialgebraic set is a finite union of basic closed semialgebraic sets @, Theorem 2.7.2]. Thus S is a 
finite union of pieces, where each piece is described by a finite number of polynomial inequalities of the type f(p) > 
with / : H' R a polynomial function. We conjecture that in the case of S there is just one of these pieces, i.e., that 
S is itself a basic closed semialgebraic set, see Definition Q] and Conjecture [4] A stronger version of this conjecture 
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is stated as Conjecture [5j It asserts that the polynomial functions / in the inequalities f(p) > defining S can be 
chosen to be invariant under the action of the unitary groups of the Hi- If the latter conjecture is true, then it is likely 
that the separability problem will be solved analytically in some low-dimensional cases where the Peres-Horodecki 
criterion fails. 

We denote by I\ the partial transposition operator on the space Hi, i = 1, . . . , n, computed in some fixed o.n. basis 
of that space. Thus, if p = pi ® • • • ® pi ® ■ ■ ■ ® p n then p Ti = pi ® ■ ■ ■ (g) pj ® • • • (g) p n . We denote by G the group 
generated by the pairwise commuting involutory operators IV The elements of G are the products T$ = Ilies ^ 
where S runs through all subsets of {1, 2, . . . , n}. We say that a state p on H has positive partial transposes (or that 
it is PPT) if p Fs > for all subsets S. If a state p is not PPT, we shall say that it is NPT. It is immediate from the 
definition of separable states that every separable state p is PPT. 

We consider often the bipartite case (n = 2). In that case we set Ha = Hi,Hb = H2, M — d\ and N = d 2 . We 
also set r = Ti. We say that a bipartite state p is a k x I state if its local ranks are k and I, i.e., rankp^ = k and 
rank /9g = / where pA = Ttb(p) and ps = Tr a(p)- Let p be a k x I state. If k = 1 or I = 1 then p is clearly separable 
(and PPT). More generally, it is well known that p is separable if it is PPT and (k - I) (I - 1) < 2 fulfil}. 

Due to the fact that the sets T> r , S r and 5^, are semialgebraic, they have well defined dimensions. We have computed 
the dimension of D r for all r = 1, . . . , d, see Theorem [7] We also found tight upper bounds for Dim S' r . When these 
bounds do not exceed d 2 — 1, the dimension of the ambient affine space H' , they are saturated in most cases. For 
several low-dimensional quantum systems we have determined the Dim S' r exactly, see Table HI 

As a simple consequence of these bounds, we prove that if (M — 2)(N — 2) > 1 then there exist separable states 
on M <§> N of length exceeding MN. (E.g. in 3 <£> 4, there exist separable states of lengths 13 and 14.) This question 
was raised about 12 years ago Q and remained open until now. In Conjecture I1Q[ we propose a new candidate for 
the maximum length of separable states. (A counter-example to this conjecture was recently found by K-C. Ha and 
S-H. Kye H .) Apparently, it is much harder to compute Dim S r . We were able to do that in 2 g) 2 and 2 (£> 3 systems, 
see Proposition QTJ 

We investigate the set S ro of separable density matrices p £ S all of whose entries are real. We denote by S 
the set of G-invariant matrices p £ S. If a Hcrmitian operator p is G-invariant, then p T = p and so p* = p. i.e., p 
is real. Thus, we have S G C S 1C . In Proposition ITUl we show that the set S consists of all states p which admit 
a representation as a convex linear combination of normalized real pure product states. We denote by L (p) the 
minimal number of real pure product states in such a convex linear combination. In Theorem[T3]we prove that, in the 
2 (£> N system, we always have L (p) = rank p. The dimensions of the sets S re and S G are computed in Proposition 
ITSl In analogy to the sets S' r , we introduce the sets S G = {p £ S : L G (p) < r}. For bipartite systems of dimension 
d < 16, we give in Table [TT] the lower bounds for dimensions of the sets S G for r = 1, . . . , d. In most cases it is shown 
that the bound is actually equal to Dim S G . We consider the following three separability problems (Si), (S2) and 
(53): for a given state p, decide whether p belongs to S, S IC or S G , respectively. Since S 2 S TC D S G , (S2) is a special 
case of (Si), and (S3) a special case of (S2). We prove in Proposition (TH] that (Si), (S2) and (S3) are equivalent to 
each other. 

The content of the paper is as follows. In section [TT] we prove that the sets V r , S r and S' r are semialgebraic and 
propose the two separability conjectures. In section UTTI wc compute Dim2? r , find a tight upper bound for DimS^,, 
and in some cases prove that they are equal. In section IIVI we characterize the set S and compute its dimension. 
Finally, in section [V] wc highlitc some of our results and discuss the prospects of solving analytically the separability 
problem in some quantum systems of low dimension. 

II. SOME SEMIALGEBRAIC SETS AND SEPARABILITY CONJECTURES 

The sets T>,S,T> r ,S r and S' r are compact subsets of H', and the first two are also convex. It is not hard to prove 
that all of them are also semialgebraic subsets of H 1 . 

Let us recall and show that V is a semialgebraic subset of H' . Indeed, if f(t) = X)i=o( — ^-) %c it d ~ % i ( c o = 1), is the 
characteristic polynomial of p £ H' , then p £ T> if and only if each a > 0. Recall that Cj is the sum of all principal 
minors of order i of p, and so Cj is a polynomial function on H' . Thus, T> is semialgebraic. In fact this shows that T> 
is a basic closed semialgebraic set according to the following definition. 

Definition 1 ( see J3|, Definition 2.7.1]) A subset X of a Euclidean space E is a basic closed semialgebraic set if there 
exist finitely many real polynomial functions fi : E R, i = 1, . . . , k, such that X = {x £ E : fi(x) >0, i = 1, . . . , k}. 

(We shall use this definition when E is either H or H' .) 

Since H r is a real algebraic set and T> r = V n H r , it follows that T> r is semialgebraic. We shall now prove that this 
is also true for the sets S^. and S r . 
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Proposition 2 Each set S' r is a semialgebraic subset of H' . 

Proof. Let Xi be the unit sphere in Hi and let the subset A C R r be defined by the equality Ai + • • • + A r = 1 
and the inequalities < Ai < 1 for i £ {1, . . . , r}. The spheres Xi are real algebraic sets and so is their product 
X = Xi x ■ ■ • x X n . As the set A is semialgebraic, the product 

Z = A x X r (1) 

of A and r copies of X is also semialgebraic. We shall now define a map / : Z — >• S' r . We shall write an arbitrary 
point z £ Z as z = (Aj, . . . , A r , X\, . . . , x r ), where each Xi = (xn, . . . , x in ) £ X. The function / is defined by setting 
f( z ) = Ei=i ^i\ x i)( x i\ G S' r i where |xi)(xi| = IxjiX^iil <g • • • <£> l^mX^inl is a pure product state. The assertion follows 
from 0, Theorem 2.8.8] since / is a real polynomial map and /(if) = S' r . □ 

Corollary 3 The set S of all normalized separable states on H is semialgebraic, and so are its subsets S r . 

Proof. Since S r = Sf]T> r , the second assertion follows from the first. We shall give two proofs for the first assertion. 

First proof: By a result of P. Horodecki [l2|, Theorem 1], which easily extends to the multipartite case, we have 
S = S' r for r = d 2 and the assertion follows from Proposition [2J 

Second proof: By Proposition [2 Si is semialgebraic. We can now apply the known fact that the convex hull of a 
semialgebraic set is also semialgebraic. As S is the convex hull of Si, we deduce that S is semialgebraic. □ 

Our first conjecture says that S is a very nice semialgebraic subsets of H' . 

Conjecture 4 The set S is a basic closed semialgebraic subset of H' . 

Our second conjecture is a stronger version of this one. It asserts that the polynomial functions, which occur in 
the representation of S as a basic closed semialgebraic set, can be chosen in a very special way. In this place, it is 
convenient to work with the non-normalized states. If X C H ' is a nonempty subset, we define the cone over X to 
be the subset KX = {tp : t > 0, p £ X} of H . Note that the vertex of KX, namely the origin of H , belongs to 
KX. Consequently, if X is closed and compact, then KX is closed. Note also that if X is convex or semialgebraic 
then KX has the same property, and conversely We are in particular interested in the cone KS consisting of all 
non- normalized separable states (plus the origin). 

The direct product of the general linear groups GL := YYi=i (C) acts naturally on TL via local invertiblc 
transformations, and also acts on H. Explicitly, if V £ GL and p € H then the latter action is given by (V, p) — > VpV^ . 
The local unitary group, i.e., the subgroup U = JlLi of G, also acts on the same spaces. We say that a polynomial 

function / : H — > R is invariant, if f(VpV>) = f(p) for all V £ U and p £ H. The Conjecture 2] implies that the cone 
KS is a basic closed semialgebraic subset of H . The following conjecture is much stronger. 

Conjecture 5 There exist finitely many homogeneous invariant polynomial functions fi : H R, i = 1, . . . , k, such 
that KS = {p £ H : f t {p) > 0, i = 1, ...,&}. 

Note that in the bipartite case both conjectures are true if (M — 1)(N — 1) < 2. In spite of the fact that this has 
been known for long time (albeit not stated in this way), all other bipartite cases still remain unsolved. If Conjecture 
[5] is true, then it should be possible to find analytic criteria of separability in some additional low-dimensional cases, 
say for 2® 4, 3® 3 and 2(g) 2(g) 2 quantum systems. To realize this objective, it is first of all necessary to find a practical 
method for computing the homogeneous invariants of small degree. Since U is compact, the algebra of polynomial 
invariants is finitely generated and the generators can be chosen to be homogeneous. We pose the following problem 
which would provide a simple method for computing all homogeneous invariants. 

Problem 6 Find a minimal set of homogeneous generators for the algebra of invariant polynomials f : H — > R for 
the quantum systems 2 <£> 4, 3 <g> 3 and 2 ® 2 g§ 2. 

III. DIMENSION COMPUTATIONS FOR SOME QUANTUM SYSTEMS 

Since all of the sets T> r , S' r and S r are semialgebraic, they have a well defined notion of dimension. We shall compute 
some of these dimensions in certain cases. We start with the sets T> r in which case we can ignore the tensor product 
structure of H. 

Theorem 7 We have Dim2? r = r(2d — r) — 1 for r = 1, 2, . . . , d. 
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Proof. Any p £ T> r has a spectral decomposition 

r 

P = j2^iM\, E A * = 1 ' A ^ ' ii^ii = i- (2) 

i=l 

We fix an o.n. basis \ak), k = 1, . . . , d of W and denote by U(<i) the global unitary group of H. with respect to this 
basis. This group acts naturally on % as well as on H and H' and its subset V. For U £ V(d) and p £ H, we shall use 
the notation U ■ p = UpU^ . Denote by X the set of all states p given by © with = \ai), i — 1, . . . , r. Clearly, we 
have D r = U(d) • X. Let Y C X consist of all p which also satisfy the inequalities Ai > A2 > • • ■ > A r . The stabilizer 
in \J(d) of any state p £ Y is the subgroup U(l) r © XJ(d — r) of dimension r + (d — r) 2 . Hence the dimension of the 
orbit V(d) ■ p is equal to r[2d — r) — r. Since Dim!" = r — 1, we infer that DimU(<i) • Y = r(2d — r) — 1. As the closure 
of V(d) ■ Y contains X, it must be equal to V(d) ■ X = T> r . Since U(d) • Y is a semialgebraic set, it follows that also 
Dim:D r = r(2d- r) - 1. □ 
Note that for r = d we recover the well known fact that Dim T> = Dim H' = d 2 — 1 . 



Theorem 8 For all positive integers r we have 

Dim<S;<r(l + 2^(di-l))-l. (3) 



Consequently, there exist separable states of length 



I := 



1 + 2E 4 (*-1) 



(4) 



where [Y| denotes the smallest integer > t. 



Proof. The space Z (see Eq. (p}) and the map / : Z — > S' r provide a parametrization of the set S' r . This 
paramctrization is redundant in the sense that the overall phases of the unit vectors \cti) and are irrelevant. To 
obtain a more economical parametrization of S' rl we replace each sphere AQ by its (2di — 2)-dimensional subsphere, 
X°, consisting of all unit vectors whose first component is real. Let X° = A° x • • • x Al° and Zq = A x (X°) r , and 
let fo : Zq —> S' r be the restriction of the map / used in the proof of Proposition [2j Since we still have fo(Zo) = S' r , 
it follows that DimS; < DimZ = r(l +2£\(dj - 1)) - 1. 

If t = I — 1 then r(l + 2 — 1)) < d 2 , and so Dim6>^. < d 2 — 1 = Dim6>. Thus S' r is a proper subset of S, and 

so there exist p £ S with L(p) > I. □ 

About 12 years ago the authors of raised the question whether the separable states on M®N satisfy the inequality 
L(p) _■ MN. It follows from the theorem that the answer to this question is negative. Indeed, if (M — 2)(N — 2) > 1 
then Dim S = M 2 N 2 - 1 > MN(2M + 2N - 3) - 1 > DimS' MN and so there must exist separable states on Ti of 
length bigger than MN. 

It is now easy to answer the analogous question in the multipartite case. 



Corollary 9 Assume that d\ > e?2 d n > 2 and n > 2. If L(p) < d for all separable states p on H, then n = 2 

and either d\ = di = 3 or c?2 = 2. 

Proof. By the theorem we have / < d and Eq. Q implies that d < 1 + 2 J^(<ij — 1). Suppose that n > 2. Then 
the function f(d±, • • • , <i„) = dj — 1 — 2 X^:(<i — 1) i s strictly increasing as a function of a single variable As 
/(2, . . . , 2) = 2" — 2n — 1 > 0, we have a contradiction. We conclude that n — 2, and the corollary follows easily. □ 

On the bipartite system M ® N, if M = 1 or N = 1, all states are separable and so £y = <S r = 2? r and DimP r = 
r(2A^ — r) — 1 by Theorem [7J In Table HI we give the dimensions of the sets S' r for several small values of n and the 
di- As iS{ C 5^ C 5g C • • • , we have Dim S[ < Dim S' 2 < Dim S' s < ■ ■ ■ . 

Let us sketch the proof of the results stated in Table H By Theorem [5] we have Dim KS' r < (2M + 2N - 3)r. If r 
is in the initial range (as secified in Table |T]), we claim that the equality holds. Define the map <p : Ha x ~Hb — > H by 
ip(a, b) = \a)(a\ <E> \b)(b\ and observe that the rank of its Jacobian matrix does not exceed 2M + 2N — 3. Next define 
the map $ r : (Ha x "Hb) 1 ' H by setting 

r 

$ r (ai, 61, . . . , a r , 6 r ) = 2J ^(ffli) h)- (5) 
i=i 

The image of $ r is exactly the cone KS' r . Since $ r is a smooth map, the dimension of K S' r must be greater than or 
equal to the maximum rank of the Jacobian matrix, J[$ r ], of the map $ r . Hence, in order to prove the claim it suffices 
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to find a point p r := (ai, 61, ... , a r , b r ) £ (Ha x "Hs) r such that the rank of J[Q r ] at p r is equal to the upper bound 
(2M+2N-3)r. Note that for r > 1 we have J[$ r ] = J[$ r _i07r]+J[^07r / ], where tt : {U A x1-L B ) r -> (W^ xH B ) r_1 and 
7r' : ("Ha x Hb) t ~~ ^ Ha x %b are the projection maps sending the point (ai, 61, ... , a r , fo r ) to (01, 61, ... , a r _x, &r-i) 
and (a r , 6 r ), respectively. Consequently, we have J[$ r ] = J[$ r _i o tt] + J[<f> o tt'] and so 



rankJ[$ r ] -rankJ[$ r _i o tt] < rankJ[0o7r'] < 2M + 2N - 3. 



(6) 



Thus, it suffices to prove the above claim only for the maximal value, r max , of r in the initial range. We have done 
that numerically for all cases in the table by choosing a random point p r and evaluating the rank of J[$ r ] at p r when 
t = r max . We used the same method to prove that Dim 5/ = M 2 N 2 — 1. The multipartite cases (those with n > 2) 
were treated similarly. 

When M = 2 and N > 1 we have I = 2N and the case r = I — 1 = 2N — 1 is exceptional. The inequality from 
Theorem E] tells us that Dim KS 2N _ 1 < 4iV 2 — 1. We claim that the stronger inequality Dim KS! 2N _ 1 < AN 2 — 2 is 
valid. Indeed, if p £ KS' 2N _ 1 then L(p T ) = L(p) < 2N—1 and so we must have det p = as well as det p r = 0. Hence, 
the set KS' 2N _i is contained in each of the two irreducible hypersurfaces det p = and det p r = 0. Consequently, its 
codimcnsion in H must be at least two. This proves our claim. Note also that if r = I = 2iV then the trivial inequality 
Dim KS' 2N < AN 2 is stronger than the one provided by Theorem [8] In these two cases we again use the ranks of the 
Jacobian matrices to prove that these improved bounds are attained. 

The use of randomness can be avoided with additional effort. In the bipartite case with M = 2 and 2 < N < 9 we 
set 



|a fc ) = |0) + (fc-l)|l>, 

k) = |0) + (fc-iV0i|i>, 

|a 2 jv-i) = |0) + {N- 

M - |o), 



\b k ) = |0) + |fc-l), 
= |0) + \2N - k 
\b 2N -i) =i|l> + |JV- 
|6ajv> = IJV — 1>, 



l) + \2N-k), 



1,. 

N- 



,2N-2- 



where i is the imaginary unit. When iV = 2we can use the same expressions except that we set I&4) = |0) + |1). Then 
one can verify that J[$ r ] at the point p r := (ai, 61, ... , a r , b r ) has rank (2A^ + 1)?' for r < 2N — 1, and has ranks 
4N 2 — 2 and 4N 2 when r is 2N — 1 and 2N, respectively. Hence, Dim K S' r must be equal to (2JV+ l)r for r < 2N — 1, 
it is AN 2 - 2 for r = 2N - 1, and AN 2 for r > 2N. 

As another example, we consider the bipartite case M = 3, N = 4 in which case r max = 13. We shall select an explicit 
point P13 where the rank of the Jacobian matrix J[$i3] is equal to the upper bound (2M + 2N— 3)r max = 11.13 = 143. 

Let us define the following 14 vector pairs (|afc), \bk)) in 3 ® 4: 



l a l) 


= |0) 






61) 


= |0) 




|a 2 ) 


= |0> + 


11) 




62) 


= |0) + |1) 




«3) 


= |0>- 


11) 




63) 


= 10) + |2> 




a 4 ) 


= |0> + 


ill) 




&4> 


= |0) - |3> 




a 5 ) 


= |0) + 


|2> 




65) 


= |0) + |3) 




a 6 ) 


= |o>- 


|2) 




be) 


= |0) + |1) + 


13) 




= |0> + 


|1H 


-|2> 


h) 


= 10) - |2) 




a 8 ) 


= |0>- 


|l>-t 


-|2) 


h) 


= |0)-<|1) 




a 9) 


= |0> + 


(1 + 


i)\2) 


bs) 


= |0)+*|1) 




OlO, 


= |0)- 


H|l) 


-|2) 


ho 


= |o> + |i>- 


H2> 


an 


= |0)- 


Hi) 


+ i|2) 


611 


= |o> + |i>- 


H\2) 


012 


= |o)- 


H|l) 


+ |2> 


fol2 


= |0)+i|2> 


+ 13) 


|ai3 


= |0)- 


H|l> 


+ i|2) 


bia 


= 10) + |2) - 


H3> 


ai4 


= |o>- 


-ill) 




fol4 


= |o)-|i). 





Set p r := (01, a r ,b r ) for r = 1, . . . , 14. We have verified that the rank of ^[$13] at the point P13 is 143, and 

the rank of J[ < I ) i4] at the point P14 is 144. Since Dimi^^ = 132, each neighborhood of the point $13(^13) must 
contain infinitely many separable states of length 13. A similar property is shared by the point <f>i4(pi4). However, it 
remains unclear whether the state $13(^13) has length 13, and the state $14(^44) the length 14. 

It is tempting to conjecture that the equality Dim S' r = DimS' r+1 implies that S' r = S' r+1 . However, it was shown 
very recently that when M = N = 3 there exist separable states of length 10. Hence, in this case Sg C S' 10 while 
Dim<Sg = Dim<S{ = 80. In view of this example, we shall propose a conjecture only for bipartite systems with M = 2 
and N arbitrary, in which case I — 2N . 



Conjecture 10 LetTi = T-La®T~Lb be the Hilbert space of a bipartite quantum system. //Dim 7^ = 2 and Dim 
N , then for any separable state p onH we have L(p) < 2N . 
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TABLE I: The dimensions of the sets iSv of separable states of length at most r in di ® • • • ® d n system. The dimension increases 
with r and reaches the maximum value d 2 — 1 for r = I, the bound defined by Eq. Q, and remains constant afterwards. 



d, (7 


niiLiai i aiiiic 


hVvcicii-i+i rvn a 1 no co 
_L!jA_Ct;JJ LlLHleLl L,di3t; 


d 2 — 1; r > I 


9 1 


9- r — 1 




^- r > 9 


9 TV- (~\ <~ N <■ CA 


f9!V _l IV — 1 ■ r <r 'IN — 1 

IZiV ± 1/ _L , / \ £il V ± 


AN 2 — "\- r — 9W - 1 


4/V 2 — 1 ■ r > 9/V 

V _L , / ^ Zi V 


o,o 


q„ _ i . r . q 




an- r > Q 


3 4 


llr — 1' r < 14 




143' r > 14 


3,5 


13r - 1; r < 18 




224; r > 18 


4,4 


13r - 1; r < 20 




255; r > 20 


2,2,2 


7r - 1; r < 10 




63; r > 10 


2,2,3 


9r - 1; r < 16 




143; r > 16 


2,2,4 


llr - 1; r < 24 




255; r > 24 


2,2,2,2 


9r - 1; r < 29 




255; r > 29 



By Theorem [U the conjectural bound 2N is the best possible. So far, it was known that this conjecture is true in the 
case of two qubits, but all other cases were open. We have proved recently [f| that it is also true for the qubit-qutrit 
system. We note that, since a separable state of rank r < 2 has length r, we have Si = S[ and 52 = S' 2 - On the 
other hand recall that Sa = S contains an open ball of H' centered at the normalized identity matrix {l/d)Id, see 
[7Jj]. Consequently, DimS^ = d 2 — 1. The dimensions of the sets S r for 2 < r < d are not known. We shall compute 
them in the two smallest bipartite cases 2 <g> 2 and 2 CS> 3. The computational method from the proof below can be 
used to compute the dimensions of the sets V r consisting of all bipartite PPT states of rank at most r. (These are 
semialgebraic sets because V r ~T> r C\ T> T .) Note that in the special cases that we consider here, we have iS r = V r due 
to the Peres-Horodecki criterion. 



Proposition 11 In 2 ® 2 we have Dim S3 = 14. In2®3 



Dim<S r = 26, 31, 34 for r = 3, 4, 5, respectively. 



Proof. It is convenient to work with non-normalized states and so we shall use the cones KS r and KT> r instead of 
S r and T> r . Any non- normalized state p of rank at most r can be written as p — C^C, where C is an r x 2N matrix 
whose diagonal entries are real and those below the diagonal are 0. Let us denote by X the real vector space of such 
matrices C. Note that Dim AT = r(4A^ — r). Thus we have a surjective map g : X — > KT> r defined by g(C) = C'C. 
We select a point Co € I such that the state po := g(Co) € V r and, subject to this condition, the rank r of the 
Jacobian matrix Jo of g at Co is maximal. By continuity, there is an e > such that the rank of the Jacobian matrix 
is equal to r and g{C) G V r for all C G X in the open ball ||C — Co II < £• The image of this ball is a submanifold 
of H' of dimension ro. Since this image is contained in KV r = KS r , we conclude that Dim KS r > ro. Hence, if 
ro = r(4A^ — r) = DimA^P r we can deduce that Dim6> r ~ DimP r . This is indeed true in the cases when (N,r) is 
(2, 3), (3, 3), (3, 4) or (3, 5). The matrices Cq for these four cases can be chosen as follows: 



2 
110 
11 



[la h], 



1 

10 10 

1 

1 1 



1 

10 10 

1 

1 2 

1 1 



(7) 



We conclude this section with an interesting example. 



□ 



Example 12 Consider the separable M x N state a = J2iLi l*)(*l ® Pi with the p; > 0. Evidently, L(<r) = ranker = 
MN. If |a, b) is a product vector and p = a + \a, b)(a, b\, then rankp = MN and we claim that L(p) = MN. There is 
a unique pi > such that Oi := pi — p.;|6)(6| > and ranker,; = A^ — 1. Then 



M 



M 



shows that L(p) < MN, and the claim follows. 




(8) 

□ 
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An interesting problem is to characterize the separable states a, say on M <S) N , such that, for every product vector 
| a, 6) on M <g> N, the state p = a + \a, b)(a, b\ satisfies the inequality L(p) < L(a). The cases where a has the maximal 
possible length are of course trivial. The above example shows that there exist a with the above stated property which 
are not of this trivial type. Since L(p) > rankp, we infer that in Example [T2l we always have L(p) = L(a) = MN. 

We claim that there exist p and a as above such that L(p) < L(a). First, we show that the sum of two separable 
states may have smaller length than one of the two summands. Assume that (M — 2)(N — 2) > 1. Then by Theorem 
[5]wc can choose a separable state a such that L{a) > MN. We can choose large t > such that p = tl — a > is 
separable. Then L(p + a) = MN < L(<r). Second, let k = L(p) and choose a decomposition p = ^2 i=1 |ai, bi)(ai, bi\. 
Let o~ r — a + |oi, &i| + • • • + \a r , b r ){a ri b r \ for r = 0, 1, . . . , k. Since L(ofc) < L(o~), there exists an r (0 < r < k) 

such that L(o~ r ) < L(cr r _i). This proves our claim. For analytic examples of such states in 3 (g> 3 see the very recent 
preprint (Toj . 

However, it is not known whether in 2 (g) N there exists a separable state a and a state p = a + |a, b)(a, b\ such that 
L(p) < L(a). If such states exist, then 

L(a) > L(p) > max(rankp, rankp r ) > max(rank cr, ranker 1 "). (9) 

One can show that there exists a separable state a whose length is larger than max(rank<7, ranker r ). Actually, analytic 
examples of such 3x3 states can be found in Q. However, no such example is known in 2 (8 JV. If such states do not 
exist, then it would follow that Conjecture [TU] is true. 



IV. REAL SEPARABLE STATES 



Many important states used in quantum information are real i.e., all entries of their density matrices are real. For 
this to make sense, we have to assume that we have fixed an o.n. basis in each of the Hilbert spaces Hi, i = 1, . . . , n. 
In this section we shall study the real separable states on H. We have defined in the Introduction the partial 
transposition operators i = 1, . . . ,n., as well as Tg for S C {1, . . . , n}, and the group G. Recall that H is the space 
of all Hcrmitian operators on H, and H 1 its affinc subspacc consisting of operators of trace 1. 

We say that an operator p £ H is G-invariant, and we write p G = p, if each element of G fixes p. (This is the case 
if and only if each Ti fixes p.) We introduce the following real spaces 

H lc = {peH:p*=p}, (10) 
H G = { P eH:p G = p}, (11) 
H' G = {peH':p G = p}. (12) 

Thus H le is the space of all real symmetric operators on H. It is easy to see that H C H re and that H G can be 
identified with the tensor product over R of the spaces of real symmetric operators onHj, t = l,...,n. In particular, 
it follows that 

Dimff G = n( di 2 +1 ). (13) 

For convenience, we also set 

2? rc = £>nir°, s ie = snH ie , s g = sdh g . (14) 

As a subset of S, the set S G is defined by the equations p Ti = p, i = 1, . . . , n. 

Let \di) £ Hi, i = 1, . . . , n, be nonzero real vectors and \<p) = \ai, . . . , a n ) the corresponding real product vector. 
Then we say that \<fi)(4>\ is a real pure product state. These states are obviously G-invariant, i.e., |</))(0| rs = \4>)(4>\ for 
all S. We say that a state p is separable over R if it belongs to the convex hull of the set of real pure product states. 
Consequently, if a separable state p is separable over R, then we must have p G = p. 

For instance, for the two-qubit real separable state p = |00)(00| + |11)(11| + |^X# with \ip) = (|01) + |10>)/V2, we 
have p T 7^ p. Therefore p is not separable over R. For an explicit expression of p as the sum of four (complex) pure 
product states see 0, Eq. (137)]. 

Our first result is that the above necessary condition is also sufficient. 

Proposition 13 Let p be a separable state and I = L(p). If p G = p then p is a sum of 2 n l real pure product states, 
and so p is separable over R. Consequently, S G is the set of all states which are separable over R. 
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Proof. We have p = £) fc=1 \tfj k )(ip k \, where |V>fc) = \a k i, ■ ■ ■ , a fc „). We define 

\b kj ) - (M + K))/V2, \c kj ) = i{\a kj ) - \a* kj ))/V2. (15) 
It is easy to verify that \a kj )(a kj \ + \a* kj )(a* kj \ = \b kj )(b kj \ + \c kj )(c kj \. Since p Ts = p for all S C {1, . . . , n}, we have 



2"/0 = ^ ^ |fflfci, • ■ ■ ,a kn )(a k i, a kn \ r 
s fe=l 

= ( \ a ki)(aki\ + |afei>( a feil J ® • • • ® ( \a kn )(a kn \ + \a* kn ){a* kn 



k=l 
I 



= (\b k l){b kl \ + \c kl ){c kl \^ «•••(» ^|6fc„X*fcn| + |CfcnXcfcnl)' ( 16 ) 

Since and Iciy) are real, this completes the proof. □ 
For p G KS G wc denote by L G (p) the smallest integer fc such that p can be written as a sum of A: real pure 

product states. We also say that L G (p) is the length of p over R. Finally, for any positive integer r, we set 

S G = {p £ S G : L G (p) < r}. It is immediate from the definitions that the sets £> ro , S G and S G are semialgebraic. 
In general, for p <E <S G we have L(p) < L G (p). It is an open question whether the equality always holds. E.g., we 

do not know whether p in Proposition [13] is a sum of I real pure product states. The next theorem shows that the 

equality L(p) = L G (fi) holds in any 2 g) N system. 

It was shown in Theorem 2] that a 2 x ]V PPT state p with p T = p is separable. Moreover, such p admits 

a decomposition p = Yli=i \ a i> ^i)( a ii °i\ with r = rank/3 and all |aj) real. For real p we have the following stronger 

version of this result. 

Theorem 14 In any 2 ® N system, every G-invariant state p is separable over R. Moreover, for such p we have 
L G (p) = rank p. 

Proof. To prove the first assertion, it suffices to consider the case where p is a 2 x N state. By the theorem cited 
above, p is separable. Thus, p £ S and the assertion follows from Proposition [131 

To prove the second assertion, we again may assume that p is a 2 x N state. It follows that r := rankp > N. Since 
p T = p, we have 



A B 
B C 



(17) 



where A, B, C are real symmetric matrices of order N. Moreover, we can assume that C is invertiblc (see e.g. [3, 
Example 2]). By performing an invertiblc real local operation on T-Lb, we can assume that C = In- By performing 
yet another local operation onT-Ls, this time with a real orthogonal matrix, we may also assume that B is a diagonal 
matrix, say B = diag(£>o, • • • j &zv— i)- Then we have 



N-l 



p= |^)(^| + |0X0|®A', (18) 

i=0 

where \4n) = (6j|0) + |1)) <E) \i) and A' = A — B 2 . Since p > 0, we must have A' > 0. Moreover, it is clear that 
r = N + r' where r 1 = rank A'. Since A' is a sum of r' positive semidefinite matrices of rank one, and the \<pi) are real 
product vectors, the assertion is proved. □ 
The first assertion of this theorem may fail for some other quantum systems. To construct a 3 x 3 counter-example, we 
start with an UPB {|a,;, 6 2 ) : i = 1, . . . , 5} consisting of real product vectors. Then the state p = I — ^2i=i \ a i> bi)(a,i, bi\ 
is G-invariant and entangled @. One can construct similarly a 2 x 2 x 2 counter-example by using the UPB from 0, 
Eq. (22)]. 

We can compute the dimensions of S lc and S G . 



Proposition 15 We have DimS rc = - 1 and DimS G = C^ 1 ) 
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TABLE II: Lower bounds for the dimensions of the sets S G (r = 1, 2, . . . , d) of G- invariant real separable states p with L G (p) < r 
in M ® N systems with d = MN < 16. 



Ivl , IV \l 


1 9 ^ 4 ^ 7 S Q Ifl 1 1 1 9 1 ^ 1 zt 1 ^ 1 fi 

1 Z O 4 O U I O y -LU 1 -L 1Z lO 14; lO -LU 


9 1 


1 9 
1 z 


9 9 


9^78 

Z O 1 O 


9 3 


■3 711 14 If! 17 

O 111 14 1U 1 1 


9 A 


A Q 14 1Q 9^ 9fi 98 9Q 
4 y 14 ly zo zu zo zy 


9 ^ 


^11 17 9^ 9Q ^4 41 A'i 44 

O 11 1 1 ZO Zy 04 OO 41 40 44 


9 


R 1 9(1 97 "^4 41 47 ^9 ^fi fil R9 

U lO ZU Z/ 04 41 4( JZ OU Oy Ul UZ 


9 7 


7 1 ^ 9" 1 ! "^I 47 ^ R9 R8 7^ 77 8(1 89 S 1 } 

1 lO zo oi oy 4/ OO UZ UO 1 O II OU OZ OO 


9 & 


S 17 9(1 W 44 ^ R9 71 7Q 8fi Q9 Q7 1(11 1 (14 1 flfi 1 (17 

O 1/ ZU OO 44 OO UZ 11 1 y OU yz ul 1U1 1U4 1UU 1U 1 


^ 1 

O, 1 


9 /I ^ 
Z 4 


^ 9 


"3 711 14 1 £ 1 7 
(11 14 1U 1 ( 




A Q U 1Q 94 9Q W ^4 

4 y 14 ly Z4 zy oz 04 oo 


3, 4 


5 11 17 23 29 35 41 47 53 56 58 59 


3, 5 


6 13 20 27 34 41 48 55 62 69 76 83 86 88 89 


4, 1 


3 6 8 9 


4, 2 


4 9 14 19 23 26 28 29 


4, 3 


5 11 17 23 29 35 41 47 53 56 58 59 


4, 4 


6 13 20 27 34 41 48 55 62 69 76 83 90 96 98 99 



Proof. Recall that there is an open ball, say B, in H' centered at the state po := (l/d)Id such that B C S. Hence, 
the first formula follows from the facts that p Q G H' IC and Dim# /re = { d + x ) - 1, where H' Ie = H' n H rc . The second 
formula follows from Eq. (|13|) by a similar argument. □ 
In Table II we exhibit the lower bounds for the dimensions of the sets S G for several bipartite systems of small 
dimension (d < 16). In most cases we have proved that these bounds are equal to Dim S G . To be precise, we know 
that the last number (for r = d) of each item is correct since it is equal to the dimension of S G as given by Proposition 
IT5l The first number (for r = 1) is also correct, it is equal to M + N — 2. The numbers following it arc correct 
as long as the difference between the consecutive numbers is equal to M + N — 1. This follows from the fact that 
Dim S^ +1 — Dim S G < M + N — 1 for each r. (This inequality can be proved by a method similar to one we used to 
prove ©.) 

Let us state the general separability problem for arbitrary multipartite systems R = Ri <g> • • • ® R. n and its two 
special cases. 

(Si) For p G T>, decide whether p G S. 

(5*2 ) For p G V n H Ie , decide whether p G S. 

(S 3 ) For p G V n H G , decide whether p G S. 

Since H G C H re , (S3) is a special case of (5*2). Next we show that (Si) can be reduced to (5*3) at the expense of 
enlarging the dimension of the quantum system. So the three problems in (Si), (S2) and (S3) are indeed equivalent. 
Since NPT states are entangled, it suffices to consider the PPT states only. 

Proposition 16 Assume that each party Ai consists of two parties: party A± t i which has a qubit and another party 
A2,i- So, we have Ri = R\,t <g> R2,i- Let p be a PPT state on R2.1 <& ■ ■ ■ <8> T~i2,n- Then there exists a state a on the 
system H., which belongs to H G , and is such that p is separable if and only if a is separable. 

Proof. The partial transposition operator Ti is the product of the transposition operators and Y" on Hxj and 
%2,i, respectively. Similarly, for any subset S C {1, . . . , n}, we have Ts = F' s T'g. Let us define the Hermitian operator 
<7, acting on the bipartite composite system TL = Wai i,....A t „ ® T~La 2 i,...,a 2 „j 

a = J2 tY " s ® P F '^ ( 19 ) 
S 

where r = \a, . . . , a)(a, . . . , o\ai i,...,Ai nJ l a ) = (|0) +i|l))/v2) and the summation is over all subsets S of the set 
{!,..., n}. Since p is PPT, we have a> 0. If R C {1, . . . , n}, then T R T S = T RA s where HAS := (R\S)U(S\ R) is 
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the symmetric difference of R and S. Consequently, we have 



a 




(20) 



The last equality holds because when S runs through all subsets of {1, . . . ,n} so does RAS. Thus we have shown 
that a G H G . Since (o|o*) = 0, we have p = (a, . . . , a|cr|a, . . . , a). This formula and Eq. ([T9|) show that p is separable 



It is well known that the set of all normalized states, T>, has nonempty interior when viewed as a subset of the 
ambient afRnc space H' , and so Dim 2? = d 2 — 1. No such result is known for its subsets T> r = {p G T> : rankp < r}. 
First of all, we have shown that T> r are real scmialgcbraic sets and so they have a well defined dimension. Then we 
have given a simple formula for their dimensions (see Theorem [7|). 

Next consider the set, S, of normalized separable states. First we dealt with its subsets S' r = {p G S : L(p) < r}, 
where L(p) is the length of p. We showed that each S' r is semialgebraic and deduced from this fact that S itself is 
semialgebraic. We have obtained in Theorem [5] very good upper bounds for Dim S' r . These bounds are not of interest 
when they exceed d 2 — 1. However, in most (but not all) of the other cases that we have computed (sec Tabic QJ 
these bounds are saturated. A simple consequence of these bounds is the fact that there exist separable states p with 
L(p) > d := Dim'H. The dimension of the subset S r = {p G S : rankp < r}, is much harder to compute. We have 
done that for the systems 2 (g 2 and 2 ® 3 only (see Proposition [TTj) . 

We have initiated the study of real separable states. It may be surprising that such states are not necessarily 
separable over R, i.e., not necessarily expressible as a sum of real pure product states (see the example above 
Proposition [13]). On the other hand we show that, among all separable states, those which are separable over R are 
characterized by the property of being G-invariant. In addition to the standard separability problem (Si), we have 
formulated in section ITVl two variations (S2) and (53) which ask to decide whether a state p belongs to the set of 
real separable states iS re or the set of G-invariant separable states S G , respectively. We have shown that all three 
separability problems are equivalent to each other. 

Last, but not least, we have proposed a method of solving the standard separability problem (Si) in some low- 
dimensional quantum systems (see Conjectures 0] and [5]). Since these very natural conjectures are valid in the two 
cases where the separability problem has been solved, namely 2 £g> 2 and 2 (g) 3, we are hopeful that they may lead to 
eventual analytic solution of the problem in some additional cases. One possibility is to use the Jacobian matrix of 
the map $ r , sec Eq. ([3]), where r is chosen so that = S. In that case, this matrix must have deficient rank at 
the points mapped to the boundary of the cone KS. This may help us to find the polynomial equations defining the 
boundary of S. 

The separabilty problem (53) is the real analog of the standard separability problem (Si). Due to the fact that the 
set S has much lower dimension than S, it is very likely that in concrete cases it would be much easier to solve (S3) 
than (Si). As we have shown that (5*3) has a very simple solution for 2 <g> N systems (see Theorem [H)) . the smallest 
open cases are 2 (8 2 ® 2 and 3 <g> 3. 
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